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Abstract 



Let rio be a polygon in R^, or more generally a compact surface with 
piecewise smooth boundary and corners. Suppose that il^ is a family of 
surfaces with C°° boundary which converges to Qq smoothly away from 
the corners, and in a precise way at the vertices to be described in the 
paper. Fedosov [4], Kac [5] and McKean-Singer [10] recognized that 
certain heat trace coefficients, in particular the coefficient of are not 
continuous as e \ 0. We describe this anomaly using renormalized 
heat invariants of an auxiliary smooth domain Z which models the 
corner formation. The result applies both for Dirichlct and Neumann 
conditions. We also include a discussion of what one might expect in 
higher dimensions. 

1 Introduction 

Let C be a domain with smooth boundary, or more generally, any two 
dimensional compact Riemannian manifold with smooth boundary. The 
Laplace operator with Dirichlet boundary conditions has discrete spectrum 
{Aj} and corresponding eigenfunctions {(pi}- The fundamental solution for 
the Cauchy problem for the heat equation has Schwartz kernel 



this converges in C°°((0, oo) x x Q) and is smooth up to t = away from 
the diagonal of O x fi. The so-called heat trace is the function 



oo 




(1.1) 
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this has an asymptotic expansion as t \ of the form 

oo 

TtH^ (1.2) 

j=0 

Each coefficient aj is a sum of two terms: an integral over of some uni- 
versal polynomial in the Gauss curvature K of the metric and its covariant 
derivatives, and an integral over dO, of another universal polynomial in the 
geodesic curvature k, of the boundary and its derivatives. Precise formulae for 
these polynomials are extremely complicated (and mostly unknown) when 
j is large, but the first few are quite simple: 



and 




Here and elsewhere, | • | refers to either area of a domain or length of its 
boundary, as appropriate. 

Almost all of this remains true if the boundary of is piecewise smooth. 
More precisely, assume that dQ is a finite union of smooth arcs, 7^, i = 
1, . . . ,k, where (counting indices mod k) ji meets "fi+i at the vertex pi with 
an interior angle G (0, 27r). In fact, the only modification in the state- 
ments above is that the heat trace coefficients may now include contributions 
from the vertices. The formula for and oi are the same as before, but 
now 

12vr ^ 4. j 2Anaj 

The term in parentheses on the right now equals 27rx{^) — Yl^=i{'^ ~ '^j)- 
That the coefficient 02 contains an extra contribution from the vertices was 
already known to Fedosov [1] (who was studying Riesz means of the eigen- 
values on polyhedra of arbitrary dimension) and to Kac [5], although the 
precise simple expression here was obtained by Dan Ray (this is referenced 
by Kac and also later by Cheeger [2] , but apparently Ray did not publish his 
result). A particularly transparent derivation of this corner term appears in 
a paper by van den Berg and Srisatkunarajah [1]. 

The heat trace anomaly in the title of our paper is the discrepancy be- 
tween the heat coefficients in the smooth and polygonal settings. More 
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specifically, it refers to the fact that at least one heat invariant is not con- 
tinuous with respect to Lipschitz convergence of domains. To phrase this 
more precisely, let 0^ be a family of surfaces with smooth boundary which 
converge to a piecewise smoothly bounded domain as e ^ 0. We think 
of as with each corner 'rounded out' slightly, but will give a precise 
formulation in the next paragraph. Denoting the heat trace coefficients for 
Of by aj{e), it will be clear from this definition that 

lim aoie) = lim ( / dA. + / n.ds 

/ KodAo + Y Kpds + y^(7r - , 

where and are the Gauss curvatures of g^, and the geodesic curvatures 
of d^le for every e > 0, respectively. The anomaly is simply that this formula 
does not agree with the expression (|1.4p . The goal of this paper is to provide 
a simple explanation for the disagreement between these two expressions. 

We now explain the desingularization more precisely. For simplicity, 
suppose that Oq and 0^ all lie in some slightly larger ambient open surface 
0, and that the metrics on are all extended to metrics (still denoted 
g^) on this larger domain. We assume that this family of metrics converges 
smoothly on il.. Let p be a vertex of and consider the portion of 0^ 
in some ball of fixed size around p, Bc{p) H il^. Our main assumption is 
that the family of pointed spaces {Bc{p) n 0^, e"^ge,p) converges in pointed 
Gromov-HausdorfF norm, and smoothly, to a noncompact region Z C 
with smooth boundary, such that at infinity, dZ is asymptotic to a cone 
with vertex at and with opening angle a, the same angle as at the vertex p 
in {QQ,g). Note that this is actually pointed Gromov-Hausdorff convergence 
for the ambient space {Q,,g^,p). 

Note that this definition implies that the distance between p and d^l^ 
is bounded above by a constant times e, and that is a small perturba- 
tion, which decreases with e, of the rescaling of the standard flat metric on 
Z n -Bc/e- For convenience we assume in the rest of this paper that the con- 
stant c equals 1. Thus the basic assumption is the existence of a smoothly 
bounded asymptotically conic region Z in the plane such that e~^{Qef^Bi{p)) 
converges to Z. 

This definition is a very special case of a more general desingularization 
construction explored carefully in [11] and [12] for the case of degeneration 
to spaces with isolated conic singularities, and in greater generality in [8]. 
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The goal in these first two papers, as here, is to analyze the behaviour 
of the heat kernel under this degeneration process. That analysis is quite 
involved, although it yields much sharper results than can be obtained by 
the present more naive methods. However, one motivation for the present 
paper is to show how some very simple rescaling arguments, which are only 
slight generalizations of ones used (in substantially more sophisticated ways) 
by Cheeger [2], already yield some interesting results. 
Now consider the function 



which is smooth on the interior of the quadrant Q = {t > 0,eo > e > 0}; 
our main theorem concerns its precise regularity at the corner t = e = 0. 
This will be decribed in terms of its regularity on the parabolic blowup of Q 
which we denote Qq. This space is diffeomorphic to Q away from the origin, 
but has an extra 'front face' F replacing the point (0, 0) which encodes all 
the directions of approach to this point along parabolic trajectories. It is 
described more carefully in §2 below. One of the goals of this paper, in fact, 
is to advertise the utility and naturality of this blowup construction. 

Theorem 1.6. Let (0^1 <7e) be a family of smooth surfaces with Riemannian 
metrics which converge in the manner described above to a surface with 
piecewise smooth boundary (r2o;5o)- Then the function G{t,e) lifts to Qq to 
be polyhomogeneous conormal at all boundary faces and corners. 

Recall that polyhomogeneity means simply that the lift of G has asymp- 
totic expansions at all boundary faces and product type expansions at all 
corners. The existence of such expansions somehow normalizes our problem. 
Indeed, the heat trace anomaly is simply the fact that the limit as e \ 
of the second asymptotic coefficient 02 (e) in the expansion as t \ is not 
the same as the second asymptotic coefficient of the heat expansion for Oq. 
The front face F of Qq separates where these limits are taken (first t ^ Q 
then e — > vs. the other way around), and this extra space allows for the 
existence of a function which interpolates between these two values. Our 
second main result describes this function. 

Theorem 1.7. There is a function C2{t) defined along the front face of Qq, 




(1.5) 
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which is smooth in the rescaled time variable r = t/e'^, and satisfies 



iim 02 (t) = , ana 

r\0 6 



k ^2 ^2 -, k 



hm C2(r) = ^ ' + > ^ > (vr - a 

r/oo ^ ^ 6 ^ 247ra,- 127r 

/is explicit form includes the finite part of a divergent expansion: 

k k 

C2(r) = ^ + ^ f .p. [ H^. ^r, z, z) dz - ±- ^(^ - a,), 

where Zj is a noncompact region in the plane which models the collapse at 
the corner. 

Remark 1.8. When VLq is a triangle (or indeed, any polygon in the plane), 
the first and third terms in the formula for liniT-^oo ^^2(7") cancel, and we 
obtain Ray's original formula 



lim C2(t) = 02(0) = Y 



247ra,- 

This interpolating function C2(t) therefore 'explains' the heat trace 
anomaly, or alternately, the anomaly is caused by the renormalized heat 
trace on the complete space Zj. We also discuss some of the other coeffi- 
cients in the asymptotic expansions for the lift of G at the various boundary 
faces and corners of Qo- 

Finally, we note that the behaviour of spectral quantities under 'self- 
similar smoothing of corners' in two-dimensional domains has been consid- 
ered elsewhere. In particular, Dauge, Tordeux and Vial [3] have carried out 
an extensive analysis of the asymptotic behaviour of solutions of An = / on 
such a family of domains. 

This paper is organized as follows. In §2, we recall some preliminary 
facts about parabolic blowups and scaling properties of heat kernels and the 
standard parametrix construction for heat kernels. The proofs of the two 
theorems are then presented in §3. In §4 we indicate the minor modifications 
needed to prove the analogous result for Neumann boundary conditions; the 
statement of the main theorem in that setting will be given there. Finally, 
in §5 we include, at the referee's request, a brief discussion of what can be 
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done in higher dimensions. The results there are less explicit, but the proofs 
carry over fairly directly. 

The authors wish to thank Lennie Friedlander for bringing this problem 
to their attention, and for making some very useful remarks on the first 
version of this manuscript; the first author is also grateful to Gilles Carron 
and Andrew Hassell for some helpful comments. 

2 Preliminaries 

In this section we collect the requisite facts and tools: the behaviour of 
the heat kernel under scaling of the underlying space, a review of parabolic 
blowups and poly homogeneity, and a slight modification of the standard 
parametrix construction for heat kernels. 

2.1 Heat kernels and dilations 

The heat kernel transforms naturally under dilations of the domain, or equiv- 
alently, of the metric. Let (M, g) be any complete Riemannian manifold with 
smooth (or piecewise smooth) boundary, and denote by {t, z, z') the min- 
imal heat kernel for the Laplacian with Dirichlet boundary conditions on M. 
This is a smooth function on the interior of x M x M with well-known 
regularity properties at the various boundaries and corners. 

We seek to relate this heat kernel with the one for the same manifold 
M but with rescaled metric gx = \^g, A G M"*". This will be applied when 
M C M^, 5 is the induced Euclidean metric, and we relate its heat kernel to 
the one for AM, the image of M under the dilation Dx : ^ 'M? , z ^ Xz. 
The pullback of the Euclidean metric from AM to M is simply A^ g. 

Proposition 2.1. The heat kernels on M and AM are related by the formula 
H^'''{\H,Xz,Xz')X^ = H^'it,z,z'). 



Implicit in this formula, we are parametrizing points in AM with points 
in M via -Da- To prove this, observe that the heat operator dt — 
on M transforms homogeneously with respect to the parabolic dilation 
{t,z) I— > (A^i,Az). Hence, the expression on the left satisfies the heat equa- 
tion; the additional A^ is the Jacobian factor accounting for the fact that 
H'^^^ {{), w, w') = 6{w — w') is homogeneous of order —2 in two dimensions. 
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2.2 Parabolic blowup 

The parabolic dilation Dx(t,e) = (A^t, Ae) motivates the introduction of 
parabolic blowup Qq of the quadrant Q := [0,oo)t x [0,eo)e at (0,0). This 
space is defined as follows. As a set, Qq is the disjoint union of Q \ {(0,0)} 
and the orbit space F = (Q \ {(0, 0)}) / ~, where (t, e) ~ {t' , e') if {t' , e') = 
Dx{t,e) for some A > 0. More concretely, F is diffeomorphic to a closed 
quarter-circle; it is also identified with the set of all equivalence classes of 
parametrized curves 7(5) = (t(s),e(s)) with lims\o7(s) = (0,0), e{s) = 
0(t(s)'^), and where 

7 ~ 7 <;=^ hm — — / = 1. 

' ' t{s) I t{s) 

The curves t = re^ (parametrized by s 1-^ (ts'^, s)), t > 0, provide represen- 
tatives of each equivalence class except the one represented by the t axis. 
There is a unique minimal C°° structure on Qq for which the lifts of smooth 
functions from Q and the parabolic polar coordinates r = Vt + e"^, t/r"^ and 
e^/r^ are all smooth. We label the faces of Qq as follows: F is the new 
front face, and L and R are the left and right side faces (the lifts of t = 
and e = 0, respectively). There is a smooth 'blowdown' map [3 : Qq — > Q 
defined in the obvious way. 

It is usually more convenient to use projective rather than polar coordi- 
nates. There are two such systems, 

(T,e), T = t/e^, and (t,r?), ri = e/Vt, 

which are valid away from R and L, respectively. Thus, for example, r is 
an 'angular' variable which vanishes on L, and in this coordinate system 
F = {e = 0}. 

Parabolic blowups are described in detail and much greater generality in 

2.3 Poly homogeneous conormal functions 

Let M be a manifold with corners. A class of functions which is the natural 
replacement for (or at least just as good as) the class of smooth functions 
is the class of poly homogeneous conormal functions. We refer to for a 
detailed exposition, but review a few facts about these here. 

First recall the space of all smooth vector fields on M which are 
tangent to all boundaries of M. If i^i, . . . , Hk are boundary hypersurfaces 
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of M meeting at a comer of codimension fc, with boundary defining func- 
tions respectively, and local coordinates y = {yi, ■ ■ ■ ,yn~k) on 
the corner, then Vb is spanned over C°°{M) locally near this corner by 
{xidxi, • • • , Xf^dxi, , dy-^ ) • • • ) }• 

A function (or distribution) u is said to be conormal if it has stable 
regularity with respect to V;,. In other words, there exists a A;-tuple of real 
numbers . . . , fXk so that 

yi...FfnG xf ...xf L°°(M), V£ andVyjGVfo. 

(In particular, the Hi are independent of £ and the Vj.) Examples include 
monomials x^^ . . . rr^'' for sj £ C, as well as products of arbitrary powers of 
I log Xj I . (This definition is slightly inaccurate since it omits the distributions 
supported at the boundary, i.e. delta sections and their derivatives, which 
are also conormal, but suffices here.) The special subclass with which we are 
interested consists of the functions with asymptotic expansions in terms of 
powers of the boundary defining functions and nonnnegative integer powers 
of the logs of these defining functions, with coefficients which are smooth 
in all other variables. The expansions are formalized using the notion of 
an index set /. This consists of a countable sequence of pairs (a, N) £ 
C X {NU {0}} such that for each A £ R, Rea > A for ah but a finite 
number of these pairs. Now, the conormal function u has a poly homogeneous 
expansion near a corner of codimension k if there are k index sets Ii, . . . ,1^ 
so that 

"~ E E a:raogxi)^i...x^'=(logXfc)^'=a„,,(y), 

{a,,N,)eIjij<Nj 

where each coefficient function is C°°. Note that since u is already 
assumed to be conormal, this expansion may be differentiated. 

The polyhomogeneous functions on Q and Qq with which we shall be 
concerned are quite simple. None of them will have log terms in their expan- 
sions, and the exponents are (not necessarily nonnegative) integers. Thus, 
for example, near L a polyhomogeneous function u will have expansion in 
powers of t with coefficients smooth in e; near F in terms of either of the 
projective coordinate systems, it has an expansion in powers of e with coef- 
ficients smooth in T, or equivalently, in powers of t with coefficients smooth 
in ij; near the corner Lf] F it will have an expansion in powers of r and e, 
with coefficients now simply numbers. 

The final point to describe here is that if u is polyhomogeneous conormal 
on Q, then its lift P*u to Qq is also polyhomogeneous conormal, and 
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On the other hand, if w is poly homogeneous on Qq, then its pushforward to 
Q is always conormal, but rarely polyhomogeneous. 

2.4 Parametrix construction 

We conclude this section by reviewing a parametrix construction for the heat 
kernel, which is useful because it accurately captures the asymptotics of the 
true heat kernel as t \ 0. The construction here is slightly nonstandard, 
but is well suited for our calculations below. 

Let M be a complete Riemannian manifold, possibly with boundary, 
and suppose that M = Mi U M2 where Mi and M2 are two manifolds 
with boundary with Mi n M2 = S a hypersurface. If M has boundary, 
assume that S intersects dM transversely, and Mi and M2 are manifolds 
with corners of codimension two. Suppose further that Mj lies in a slightly 
larger complete manifold Mj, again possibly with boundary, such that for 
some neighbourhood U oi J^, Mj nU = M nU . 

Taking the heat kernels on each Mj as given, define 

2 

where Xj is the characteristic function of Mj in M. In the more customary 
parametrix construction, the Mj are relatively open in M, and Mi n M2 
is also open; the H^^i are pasted together using cutoff functions {V'j} and 
{ipj} with V'l + V'2 = 1, where suppipj C {ipj = 1}, and suppV'j C Mj. We 
are using sharp (discontinuous) cutoffs rather than smooth ones, however, 
so that we can identify certain asymptotic coefficients in the calculations to 
follow. 

Lemma 2.2. Let H^^{t, z,z') denote the true heat kernel on M, and set 

K{t,z) = H^^{t,z,z) - H^^{t,z,z). 
Then K{t, z) = 0{t°°) ast\0. 
Proof. Rewrite 

H''{t,z,z)=xi{z) [H^'^{t,z,z)-H^'{t,z,z)) 

+ X2{z) (H^'Ht, z, z) - H^'it, z, z)) + H^'it, z, z). 
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By assumption, Mj agrees with M in a neighbourhood of Mj, so that 

H^^i {t, z, z) — {t, z, z) = 0{t°°) on the support of Xj (remember that 
the small t expansions of these operators are local), and this proves the 
claim. □ 



3 Proofs of main theorems 

We have now assembled all the requisite facts and can proceed with the 
proofs of the main theorems. 

As in the introduction, let G{t,e) = Tri/^^ U P : Qq ^ Q is the 
blowdown map, then let G = (3*G. We need to analyze the behaviour of 
Q near each of the faces and corners of Qq, and for that we shall use the 
coordinates (r, e) introduced in §2.2. 

We shall make a simplifying assumption about the geometry in order to 
elucidate the proof. For each i, let Sa^ denote the sector in with opening 
angle Oj. Choose a smoothly bounded region Zi in the plane which coincides 
with Sen outside i?i/2(0)) a-^^d let Zf = Bii^{Q) n Zi. Then we assume that 
near each vertex pi, the restriction of the metric to Bi[pi)r\^e is isometric 
to the dilation by the factor e of the region Z^, which obviously lies in the 
unit ball. The result remains true in the generality with which it was stated 
earlier, but the proof requires a few more technical steps which are both 
standard and not particularly germane to the main ideas here. Furthermore, 
for notational convenience only, we assume that there is only a single vertex 
p and denote the corresponding smooth model region and sector by Z and 
S, respectively. 

Proof of Theorem We first construct a particular family of parametrices 
for the heat kernel on O,.. For any < e < eo, decompose 

fi, = U Vl' 

where = $7, n Bi{p), and 0' = 0^ \ (il^ n Bi{p)). Note that 9! is 
independent of e. Lemma 12.21 shows that 

z, z) = Xi{zW^{t, z, z) + X2{zW^'{t, z, z) + K{t, z), (3.1) 

where xi is the characteristic function of |z| < 1, X2 = 1 — Xii K is the 
error term from Lemma 12.21 hence 

G{t,e)= [ H'^{t,z,z)dz+ [ H^"{t, z, z) dz + [ K{t,z,z)dz. 
J\z\<i Jo.' Jn^ 
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We denote the sum on the right side by I + II + III, and analyze the hfts of 
these terms successively. 

By Proposition EH H^^(t,z,z') = e~^H^ {t/e^ , z/e, z' /e), so setting z = 
z' = ew, we see that 

f3*l= [ H^{T,w,w)dw. 

J\w\<l/e 

This will be the principal term, and we defer its analysis for the moment. 

Next, II is independent of e, and it is polyhomogeneous as t \ 0, with 
expansion given by integrating the standard heat coefficients aj{z) over this 
restricted domain. Hence its lift to Qq is clearly polyhomogeneous. 

Finally, by Lemma [2.21 III depends on e but decays rapidly in t uniformly 
in e. 

We now examine (3*1 more closely. Choose a smoothly bounded compact 
region W which agrees with Z in \w\ < 2, so that Z = {W D Bi) U (5 \ Bi). 
Using Lemma 12.21 again, write 

H^{t, z, z) = Xi{z)H'^{t, z, z) + X2{z)H'^{t, z, z) + Ki{t, z) (3.2) 

where Ki is the corresponding error term. Then 

(3*1= [ H^{T,w,w)dw + [ H^{T,w,w)dw 

J\w\<l -'l<|w|<l/<: 

+ / Ki{T,w,w)dw, 

J\w\<l/e 

which we write as I' + II' + III'. 

We first prove poly homogeneity of these terms away from the right face 
R of Qq. The term I' has an expansion as r \ and is independent of e, so 
/3*I' is certainly polyhomogeneous in this region. By Lemma 12.21 again. Ki 
decreases rapidly as r ^ 0, so this term is also polyhomogeneous there. Note 
too that by the explicit form of the error term in the proof of that lemma, 
and using the dilation properties of and again, Ki{t, z) = 0{\z\~°°) 
uniformly for r in any bounded set, so its integral over \z\ < 1/e is also 
bounded independently of e. 

To analyze the remaining term, set 

D{R) := / H^{l,w,w)dw. 
J\w\<R 

By Proposition EH II'(e,r) = D{l/e^) - D{l/y/T), so it will suffice to 
show that D has an expansion in powers of 1/-R as i? — > co. For this, we 
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appeal to a calculation by van den Berg and Srisatkunarajah [T], who prove 
that 

°<^) - ^ - VT^* . ^ . ,3.3) 

for some c > independent of R. Only the polyhomogeneous structure of 
the second term on the right is not completely obvious. For that, we may 
as well replace the upper limit of integration by 1/2 since the integral from 
1/2 to 1 decreases exponentially in R. Using the Taylor series for \J\ — 
at y = 0, we find that 

and this completes the proof of polyhomogeneity of for r in any bounded 
set. 

To finish the proof, we must analyze the behaviour of /3*I as r ^ oo. 
Switch to the coordinates t, 77, so e = r\\ft and r = ry~^. It is now more 
convenient to use the standard representation of the heat kernel in terms of 
the resolvent: 

= j e^^^RziX)dX, where Rz{X) = {-Az - Xy\ (3.4) 

Here T is a path surrounding the spectrum of Az, for example, the two 
half-lines Im A = ib(aReA + P), a,(3 > 0, joined by the half-circle |A| = /3, 
Re A < 0, traversed in the counterclockwise direction. In general it is a subtle 
matter to deduce the fact that has an expansion in powers of 1/r at 
large times since this depends on the fine structure of the resolvent near the 
threshold A = 0. However, in this case we already have sufficient information 
about the heat kernel on S that this is not hard. Choose a partition of unity 
{V'l) V'2} on Z such that ^1 = 1 in \z\ < 3/4 and '(/'2 = 1 in l-^l > 5/4, and that 
both Z CiW and Z D S contain the region Z n {3/4 < \z\ < 5/4}. Choose 
other cutoff functions ipj such that ipj = 1 on suppV'j- Let Rw and Rs 
denote the resolvents for Aw and As, with Dirichlet boundary conditions, 
and define the parametrix 

Rz{X) = ^iRw{X)^i + i^2Rs{X)^2. 

Then 

{Az - X)Rz{X) = 1+ [Az, ^i]iiH/(A)Vi + [Az,i^2]Rs{X)iP2 ■= I + E{X). 
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Since the singular supports of both Rw and Rs are on the diagonal, and 
the support of [Aj-f/'j] is disjoint from that of ipj, we see that E{X) is a 
holomorphic family of operators (for A G C \ M"*") which maps L2(Z) into 
Cq^{Z). Also, since —Az — A is invertible for A in this region, / + E{X) is 
also invertible there. We write its inverse as / + F{X), so that 

Rz{X) = Rz{X) + Rz{X)FiX). (3.5) 

The relationships (/ + ^;(A))(/ + F{X)) = (/ + F{X)){I + E{X)) = I imply 
that 

F(A) = -^(A) + ^^(A) + E{X)FiX)E{X), 

hence F{X) is also smoothing and maps L'^{Z) into Cq^{Z); the second term 
on the right in (j3.5p has the same mapping properties. 

Finally, the form of the expansion of F{X) for A near (away from the 
positive real axis) is the precisely the same as that of E{X), which in turn 
is the same as that of Rz{X), and hence finally as that of Rs{X). 

The decay of each of term as |A| ^ oo is straightforward, so we can write 



H^{ri-\ z,z)dz= f ( [ e-^/'^^RziX) dx] (z, z) dz 

J\z\<i/'r,Vt \Jr ) 

= [ ( [ e-^/i^tl;iRw{X)dx] {z,z)dz 
J\z\<2 \Jr J 

+ [ ( [ e-^/'^'i^^RsWdx) {z,z)dz 

Jl/2<\z\<l/riVt \Jt J 

+ I (I e-^/"^^ Rz{X)F{X)dx\{z,z)dz. 

J\z\<i/Ti^/t \Jr J 



'\z\<l/riVi 

The inner integrand in the first term on the right extends holomorphically 
to a neighbourhood of A = 0, so the contour can be moved to lie entirely 
in the right half-plane, which shows that this term decreases exponentially 
in l/rj. The second term is polyhomogeneous by the explicit analysis of the 
function D{R) above. The fact that the final term has an expansion follows 
from the existence of asymptotics of F[X) for A near 0. This completes the 
proof of the poly homogeneity of ^ on Qo- D 

Proof of Theorem (jl.7p This consists of examining the terms in the expansion 
of Q at the various boundary faces. 

First, at L, away from F we may use the variables (t, e), and 

oo 

G{t,e)^Y.^,{e)t-'+^/\ 

3=0 
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Near L n F, we substitute t = e^r to get 



g(r, e)^Yl a,{e)T-^+^/h-^+K (3.6) 

j=0 

The coefficients aj{e) are polyliomogeneous as e — > by Theorem 1 1.6[ 
At R, away from t = 0, 

oo 

G(t,e)~5]i?,(t)e^; 

j=0 

here Bo{t) = TtH^". Near F n R we use the coordinates t and rj = e/ \ft to 
compute 

oo 

Again, the coefficients Bj{t) are polyhomogeneous in i. 

Finally, near F, we use the coordinates (t, e), so the expansion is in 
powers of e, and by ()3.6p it is 

oo 

i=o 

We shall identify the coefficients Cq, C\ and C2. 

By our analysis of the terms I', II', III', II and III, we see that only II' and 
II contribute to the coefficients of and e~^. Substituting directly from 
the expansions of these two terms (using the McKean-Singer asymptotics 
on 0' for II and the first terms in the expansion of L'(l/e-^) for II'), and 
then using the definition of the finite part at e = of I, we have 

1 /lO'l a\ 1 (\d^\ 1 



— (/ KdA+ [ Kds] +f.p. [ H^{T,w,w)dw + 0{e). 
271" \Jn' Jan' J e=o Jz 



In other words, 

T \ Att on 
1 f\dn'\ 1 
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and 

C2{t) = — ([ KdA+ [ Kds] +f.p. [ H^{T,w,w)dw. 
12vr VJq' Jan' J £=o Jz 

This simplifies using the following observations: first, the area of a cir- 
cular sector of opening a and radius 1, i.e. \Clo n Bi\, equals a/2, so the 
coefhcient of e~^T~^ is just |J7o|/4vr; similarly, the sides of this circular sec- 
tor are straight lines, so |3r2oni?i| = 2, which means that the next coefficient 
is —\dilo\/8y/7r; finally, since go is flat inQoHBi, K = there, so using that 
the contribution from 'turning the corner' at p in the boundary integral is 
TT — a, we flnd that 

/ KdA+ nds = 2ttx{^o) - (vr - a). 
Jn' Jan' 

This means that 

C2(T)=f.p. / H^{r,w,w)dw + ]:xm-^^- (3-7) 

We conclude by calculating its behaviour for small and large r. Using 
the small r asymptotics, we see that 



H {T,w,w)dwr^ T 



|«;|<l/e 47r 



-T 



]- f Kds + Oier'/^y, 

Joz 



S/tF 127r 

hence the flnite part of this integral is equal (up to the factor 127r) to the 
integral of curvature on the boundary of Z, which is the total turning angle 
TT — a, so flnally, the limit of C2 as r ^ is x(^^o)/6, as claimed. 
Finally, we use the dilation one more time to calculate that 

/ H^{T,w,w)dw= / H^^^ {l,w,w) dw. 

J\w\<l/€ -'k|<l/ev^ 

Noting that e-^/r = and since converges to the sector 5 as r — > 

00, we can use the expansion ()3.3p to see that the finite part is indeed 
(vr^ — a^) /247ra. Therefore, in general, with an arbitrary number of vertices, 

^lun C,(r) = ^^Y. - 12^ 2:(- - ^^^^ 

j=i j=i 

in particular, if is a polygon, its Euler characteristic is 1, so the first and 
third terms cancel. 

This completes the proof. □ 
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4 Neumann boundary conditions 

We now briefly discuss the minor modifications needed to prove the ana- 
logues of Theorems 11.61 and 11.71 if Neumann conditions are used instead of 
Dirichlet conditions. 

A cursory inspection of the proof shows that the only real issue is to 
find an analogue of the van den Berg-Srisatkunarajah formula (j3.3p in this 
setting. This does not seem to appear explicitly in the literature, but for- 
tunately, a recent paper by Kokotov [6] contains the corresponding formula 
for the complete cone of angle 2a. Let denote the heat kernel on 
this cone. Then by Proposition 1 in [6], there exists c > such that for 
every R > 0, 



This formula is stated in [6] for fixed radius R and for the heat kernel at time 
t as t ^ 0, but because of the usual scaling properties, it holds equally well 
for fixed t, say t = I, and as the radius R ^ oo; indeed, the quantity on the 
left depends only on the ratio R/t^. The coefficients in this expansion have 
been written in a nonreduced form in order to emphasize the dependence 
on the angle 2a. 

We now observe that the cone C2a is the union of two copies of the sector 
Sa with the boundary rays identified. Alternately, let r be the obvious 
reflection on the cone C2a] then a region isometric to the sector Sa is a 
fundamental domain for this action and its image T{Sa) is the other half of 
the cone. In any case, using this, the formula for the Neumann heat kernel 
follows directly from ([331) and (gT]). Indeed, let L'^iC^a) = L\ ® L'^_ be 
the decomposition into functions which are even and odd with respect to 
r. If ti G H'^{C2a) n L^, then u has vanishing normal derivative at dSa, 
while if n € H^{C2a) H then u vanishes at dSa. Since the Laplacian 
commutes with r, the heat kernel has a 2-by-2 block decomposition: the 
upper left and lower right on-diagonal blocks are canonically identified with 
the Neumann and Dirichlet heat kernels of S^, and we denote these by 
and H^, respectively. Therefore, 




(4.1) 




whence 




(4.2) 
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In other words, in this asymptotic formula, only the signs of the odd powers 
of R are reversed from those in the corresponding formula for the Dirichlet 
heat kernel. 

It is now a simple matter to track through the various arguments in 
this paper to obtain that Qf^{T,e), the pullback to Qq of the trace of the 
heat kernel for the Laplacian with Neumann boundary conditions on O^, is 
polyhomogeneous and has the expansion 

r f , 1 l^ol ^ 1 I^J^oL 

— ( [ KdA+ [ Kds]+i.p.[ H^{T,w,w)dw + Oie). 

12vr \Jq, Jqq, J e=0 J\w\<l/e 

In particular, the coefficient C2(t) of e*^ is exactly the same as in the Dirichlet 
case. We leave the straightforward details to the reader. 



5 Higher dimensions and other generaUzations 

We have focused in this paper on two-dimensional domains in order to em- 
phasize the simplicity of the arguments and to take advantage of the ex- 
plicit nature of the various formute. There are various analogues of these 
results in higher dimensions, which we now describe briefly. These general- 
izations should have some interesting applications, which will be developed 
elsewhere. 

The simplest generalization is to consider a family of Riemannian met- 
rics on a compact manifold M such that {M^g^^) degenerates to a space 
(Mo,(7o) which has isolated conic singularities. We assume that this de- 
generation is modelled on the rescalings of a complete asymptotically conic 
space {Z,gz), i.e. such that suitable neighbourhoods of {M,g^) are (asymp- 
totically equivalent to) rescalings of truncations of {Z,gz)- The behaviour 
of the entire heat kernel for this type of degeneration was studied in great 
detail in [11], [12]. The analysis in those papers is much more general, but 
considerably more intricate, than what is done here, but one consequence 
of those results is the fact that the trace of the heat kernel for Ag^ lifts to 
Qq to be polyhomogeneous, and the terms in its expansions at the various 
faces can be determined explicitly. In particular, the coefficient C„/2(''") of 
e'^ at the front face of Qo is equal to the sum 

/ Qn/2dV + f.p. / H^{T,w,w)dw, 

J Mo t=0 JZe 
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where g„/2 is the standard heat invariant integrand for the metric qq. It is 
unhkely, however, that one could find an exphcit formula for the limit as 
T — > oo of this regularized trace, except in special cases. 

There should be a similar generalization of these ideas to the setting of 
resolution blowups of iterated edge spaces (or smoothly stratified spaces), 
introduced in [8]. This is likely to present a greater challenge, but one of the 
motivations for exploring the rescaling methods used here is to find way to 
circumvent the machinery of [11] , [12] if one is only interested in heat traces 
rather than more extensive information about the heat kernel. 

A special and very interesting case would be to find an analogue of Theo- 
rem [L7l for smoothings of Euclidean polyhedra in arbitrary dimension. The 
description of a family of 'self-similar' smoothings of an arbitrary polyhe- 
dron is not difficult and follows the scheme presented in [8] closely. However, 
in order to make this formula explicit, one would need an analogue of (|3.3p 
or ()4.ip for higher dimensional polyhedral sectors, which does not seem to 
be available. (Analogous results are known for other spectral invariants, 
however, see |1] and [2].) 
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